% Het uiteindelijke resultaat van dit script is per sector in het
% net een array met daarin per belastingsbus wat het jaarlijks kost aan PQ-problemen.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

'&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&'
'&&&&&&&&&&&   Opmaken van de kosten per sector op de verschillende belastingsbussen   &&&&&&&&&&'
'&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&'

gegevens;
aantal_bussen
aantal_sectoren

%Bedoeling van onderstaande lus is om sector per sector de kosten van die sector op de verschillende bussen te berekenen en op te slaan in een array.
%Hiertoe loopt een buitenste lus over alle sectoren, een binnenste lus
%loopt over alle bussen en sommeert de kosten van die sector op die bus (hiertoe loopt een derde lus over alle rijen van de problemen_bus_j matrix).

for i=1:1:aantal_sectoren
    naam_sector = sectoren(1,i);
    test_string=strcat('kosten_sector_',naam_sector);
    test_string=char(test_string);
    %if 1==exist(test_string)
    smatrix=zeros(aantal_belastingsbussen,2);
    smatrix(1:aantal_belastingsbussen,1)=belastingsbussen(:,1);
    for j=1:1:aantal_belastingsbussen
        nummer_belastingsbus = belastingsbussen(j,1);
        string_nummer_belastingsbus=int2str(nummer_belastingsbus);
        if exist(strcat('problemen_bus_',string_nummer_belastingsbus))==1
        naam_problemen_bus_j_array = strcat('problemen_bus_',string_nummer_belastingsbus);
        kmatrix = eval(naam_problemen_bus_j_array);
        [r,k]=size(kmatrix);
        kost_sector=0;
        for p=1:1:r
            kost_sector=kost_sector+kmatrix(p,7+i);
        end
        smatrix(j,2)=kost_sector;
        end
    end
    eval([test_string '= smatrix'])
    %end
end